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We investigate unitarity bounds in the most general two Higgs doublet model without a 
discrete Z 2 symmetry nor CP conservation. S-wave amplitudes for two-body elastic scatter¬ 
ings of Nambu-Goldstone bosons and physical Higgs bosons are calculated at high energies 
for all possible initial and final states (14 neutral, 8 singly-charged and 3 doubly-charged 
states). We obtain analytic formulae for the block-diagonalized scattering matrix by the 
classification of the two body scattering states using the conserved quantum numbers at 
high energies. Imposing the condition of perturbative unitarity to the eigenvalues of the 
scattering matrix, constraints on the model parameters can be obtained. We apply our 
results to constrain the mass range of the next-to-lightest Higgs state in the model. 

PACS numbers: 


I. INTRODUCTION 

The Higgs boson was discovered at LHC in 2012 ii, and its mass and coupling constants 
turned out to be consistent with the predictions in the standard model (SM) [^, However, 
although the Higgs boson was found, the structure of the Higgs sector remains unknown. There 
are possibilities of non-minimal Higgs sectors which contain additional Higgs bosons. They also 
can satisfy the current data. In fact, in large number of new physics scenarios, extended Higgs 
sectors are predicted. The structure of the Higgs sector is strongly connected to these new physics 
scenarios. Therefore, the Higgs sector is a probe of new physics. 

Although parameters of an extended Higgs sector are basically free, they can be constrained by 
imposing some compelling theoretical conditions, such as those of perturbative unitarity, vacuum 
stability and triviality. The requirement of perturbative unitarity |Q] is known to give a conservative 
bound on the parameters of a model, beyond which the perturbation calculation does not work. 
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It goes without saying that the parameters can also be constrained by taking into account various 
experimental results such as electroweak precision data observables in flavour physics 0, the 
Higgs boson signal strength ii and so on. By the combination of the bounds from the theoretical 
consistency and those from the experimental data, a parameter space in extended Higgs sectors 
can be further restricted. 

The theoretical constraints have been intensively studied in the two Higgs doublet model 
(THDM). Nice reviews for the THDMs have been presented in Refs. ||8|, Previously, the uni- 
tarity bound was mainly studied for the model with a (softly-broken) discrete Z 2 symmetry ||l^ 
with CP-conservation EMI]- The softly-broken discrete Z 2 symmetry is phenomenologically 
important in order to avoid flavor changing neutral currents (FCNCs)^, under which there are 
four types of Yukawa interactions [IT-^. Bounds from vacuum stability and triviality have also 
been studied in Refs, and in Refs. p6|-p^], respectively, in Z 2 symmetric THDMs. These 


theoretical bounds have been used to limit magnitudes of cross sections and decay rates of tree 


level processes | ‘K , , radiative corrections due to additional Higgs bosons to the Higgs cou¬ 
plings [^, ^ and one-loop induced processes . 

However, analyses in the most general THDM without the Z 2 symmetry is getting important as 
an effective description of more various new physics scenarios, such as supersymmetric SMs with 
non-holomorphic Yukawa couplings and also general models with CP-violation [^| which is 


required for successful scenario of electroweak baryogenesis |36-^. 

In this letter, we investigate unitarity bounds in the most general THDM without the Z 2 
symmetry nor CP-conservation. We calculate S-wave amplitudes for two-body elastic scatterings 
of Nambu-Goldstone (NG) bosons and physical scalar bosons at high energies for all possible initial 
and final states (14 neutral |l^, 8 singly-charged El) and 3 doubly-charged states). By choosing 
appropriate bases of the scattering states which are deduced using the conserved quantum numbers 
at high energies such as the hypercharge, weak isospin and the third component of the latter f^ , 
the scattering matrix is given as a block-diagonalized form with at most 4x4 submatrices. Thus, all 
of the eigenvalues of the scattering matrix can be easily evaluated numerically. The analytic result 
for the block diagonalized S-wave matrix is consistent with that given in Ref. By requiring 
that each of the eigenvalues is not too large to break validity of perturbation calculation, the model 
parameters, e.g., the masses of extra Higgs bosons and mixing angles, can be constrained. Our 


^ As an alternative way to avoid the tree level FCNCs, the aligned THDMs have also been known, where two 
Yukawa matrices are assumed to be proportional with other. 
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results can be useful to constrain parameter spaces whenever one evalnates physics quantities in 
the most general THDM. 

We then numerically demonstrate that the bound on the mass M 2 nd of the second lightest Higgs 
boson, whatever it is, is obtained by inputting the mass of the discovered Higgs boson h under the 
assumption of non-zero deviations in the Higgs boson coupling hVV with the weak boson (V = W 
and Z) from the SM valne. Cnrrently, the hVV coupling was measured with ~ 10% accuracy 
by the LHC Run-I experiment (see Ref. and references therein), and that is expected to be 
measured more accurately at future collider experiments. If a deviation is detected in the hVV 
coupling, we can obtain the constraint on the region of M 2 nd even without its direct discovery. 


II. MODEL SETUP 

A. Higgs potential 


The most general Higgs potential under the SU{2)l x U{1)y gauge symmetry is given by 


V = |<hip + m\\^ 2 \^ — + h.c.) 

+ “1“ -^31*1^1 Pi‘ 1*2P + A4|<h|<I>2p + — A5(<1>|$2)^ + h.C. 


+ 


A6|^>iP^> 1$2 + A7|$2p$1^2 + h.C. 


( 1 ) 


where [i = 1,2) are the isospin doublet scalar fields with hypercharge Y = 1/2. In general, 
1712 and A 1 -A 4 are real, while m 3 and A 5 -A 7 are complex. Thus, there are totally fourteen 
real parameters. If there is a symmetry in the potential, the nnmber of parameters is rednced. 
For example, when the potential is exact (softly-broken) Z 2 invariant under the transformation of 
(<I>i,<I> 2 ) —>• (^ 1 ,-^ 2 )) the m|, Ag and A 7 (Ag and A 7 ) terms are forbidden. 

By using the U(l)y invariance and rephasing the doublet fields, the vacuum expectation values 


(VEVs) of the two doublet helds can be taken to be real without loss of generality [^-^|. The 
two doublet fields are then described in terms of the component fields as 





^^(u T hj -j- izi') 


(* = 1 , 2 ), 


( 2 ) 


where vi and V 2 are the VEVs of <I>i and <I> 2 , respectively, which are satisfied v = ~ 

(v^Ci?)"^/^ ~ 246 GeV. By introdncing tan/3 = ^ 2 /^ 1 , two VEVs are described by v and tan/3 as 
the usual notation. In the following, both the VEVs are assumed to be non-zero, except for the 
case of the inert doublet model discussed in Appendix. 
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The stationary conditions of the scalar potential are given as follows 

dV 

= 0, {ipa = hi, h2, Zi, and Z 2 ), 


difa 

where the left hand side of the above eqnation for each is calcnlated by 


dV 

dhi 

dV 

dh2 

dV 

dzi 


= VCr, 


= VSf 


mf — + — (Aic^ + A 345 S^ + 3A^ tan /3) 


m\ - cos^ + A 345 C^ + cot 13 + 


= —ntan/3 


dV 

dzn 


= VSf. 


Imm| + — {Xlsi^Ci^ + AgC^ + A^s^) 


where we used the abbreviation of Sq = sin0 and Cg = cos 6. In Eqs. ®-i) , we introduced 

X X X X R , Re mo 

A 345 = A 3 + A 4 + A 5 , M = 


■5/3 C/3 


(3) 

(4) 

(5) 

( 6 ) 

(7) 


and 


Af^ = ReAfc, A^ = ImAfe, {k = 5, 6 , 7). ( 8 ) 

From the stationary conditions in Eq. (|^, we can eliminate m^, m^ and Imm^ in the Higgs 
potential. 

In order to calculate the masses for the scalar bosons, it is convenient to introduce the so-called 


Higgs basis [45| dehned as 




(9) 


where 


<h = 


G+ 

j=,{v + h[+iG<^) 


T = 


H+ 

^(/l'2 i/1'3) 


( 10 ) 


where G^ and G^ are the NG bosons which are absorbed into the longitudinal components of 
and Z by the Higgs mechanism, respectively. The physical singly-charged scalar state is denoted 
as H^. The three neutral states /i'^, /12 and /13 are not the mass eigenstates at this stage, which 
generally mix with each other. 


B. Mass spectrum 

We give the mass formulae for the Higgs bosons in the most general case. First, the squared 
mass of is calculated by 


m‘jj± = ^ (A 4 + Xf — cot P — Xf tan /3) . 


( 11 ) 
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Next, the mass term for the neutral scalar states is expressed by the 3x3 matrix as 


T Amass U 

•^neutral — 


where each of the matrix elements is given by 


^ Ml Ml Ml ^ 



X^l ^22 -^23 


/i' 

V Mfg M23 M33 y 


V ^3 J 


Ml = M^ + -y 


Xicfj + A2S^ H ^'®2/3 + 2(Afc^ + Xfsp)s2l3 

(Ai + A 2 — 2A345 )s2/3 ~ 2 A^(s4 ^ + cot/3) — 2 A^(s4 ^ — tan/3) 
Is = (^X^ + ^ cot /3 + ^ tan /3^ , 

y2 

[{X 2 Sp — AiC^ + A345C2 ^)s2^ + 2 A^( 2 c 2 ^ — l)c| + 2Af ( 2 c 2 ^ + l)s^] 

o 


/\,f 2 _ 

1 U 12 — 


Ml — — — {Xls2j3 + 2 Agc| + 2 A 7 SI) 

y2 

Xil = ("^5^2/3 ~ X^C^ + XjSp) . 


( 12 ) 


(13a) 

(13b) 

(13c) 

(13d) 

(13e) 

(13f) 


The mass eigenvalues mjj, {i = 1-3) are obtained by introducing the 3x3 orthogonal matrix R as 

RfkMiiRij = diag(mf7^, , ml^), (14) 


where niHi ^ ''^^2 < mHz i® assumed. Mass eigenstates for the neutral Higgs bosons are also 
defined using R as 




^ Hi ^ 

h '2 

= R 

H 2 

V ^3 J 


[hs) 


(15) 


where Hi is defined to be the SM-like Higgs boson with the mass of about 125 GeV, i.e., m^_^ ~ 125 
GeV. In the following, we represent Hi and its mass m^_^ by h and m^, respectively. In the GP- 
conserving limit, H 2 and H^ respectively correspond to the additional GP-even {H) and CP-odd 
(A) Higgs bosons. The mass formulae for the CP-conserving case will be discussed in the next 
subsection. We here note that R can be described by three mixing angles In our numerical 
analysis given in Sec. IV, the matrix elements of R are derived by inputting the mass matrix 
elements given in Eq. ([l^. 


By using Eqs. (0), S), (p.3bp, (p.3bD, and (p.dcj), the five parameters (Ai _4 and X§) can be 
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rewritten as follows 

\iv^ = Mh + {M 22 — M'^) tan^ /3 — 2 Mf 2 tan ^ tan^ /3 — ^Xq) tan /3, (16) 

\ 2 v‘^ = Mfi + (MI 2 - cot^ /3 + 2 Mf 2 cot /? + i (A|* cot^ /3 - 3Af) cot /3, (17) 

Asr;^ = Mn - (MI 2 + M"^) + 2 Mi 2 cot 2/3 + 2m^± - ^(A|^cot/3 + Af tan/3), (18) 

A 4 t'^ = + Mfg — 2m'jj± — ^{Xf cot /3 + Af tan /3), (19) 

X§v‘^ = — Mfg — ^(A|*cot/3 + Af tan/3). (20) 


Therefore, we can choose thirteen input parameters as v, tan/3, M^, MI 2 , M^ 2 ) -^ 33 ) ^‘h±^ 
A ^7 and Agg 7 . Instead of MI 2 and M^ 2 ) "'^c can take d, and ?ti| as inputs by 

^-^11 = (21) 

M22 = ( 22 ) 

^-^12 = {ml - (23) 

In addition, we rewrite MI 3 as 


MI 3 — rh^. (24) 

In the CP-conserving limit (Imm| —)• 0 and A 5 g 7 —)• 0), rh|, rh\ and a describe the physical 
squared masses for h, H and A, and the mixing angle for the CP-even Higgs state, respectively. 

Consequently, by using experimental values of v (= 246 GeV) and ruh (= 125 GeV), eleven 
input parameters can be chosen as follows 


{mjj, m\, m]j±, tan/3, sin(/3 - d), jAgl, IA 7 I, 6 * 5 , Oq, 0y}, 


(25) 


where ^ 5 , 6,7 are the complex phases of As^gj- Notice that mh is determined so as to keep nih = 125 
GeV. 


III. UNITARITY BOUNDS 

We calculate the S-wave amplitude matrix for the elastic BB' ^ B''B'" scatterings in the high 
energy limit, where the fields B, B', B” and B'" represent either Hi, H 2 , H^ or In 

this case, all the longitudinal components of the weak gauge boson states can be replaced by the 
corresponding NG boson states because of the equivalence theorem. Furthermore, only the scalar 
boson contact interactions contribute to the S-wave amplitude. Therefore, the calculation of the 
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S-wave amplitude matrix is quite simply done just by extracting the coefficient of the scalar boson 
quartic terms. 

There are 14 neutral, 8 singly-charged and 3 doubly-charged 2 body scalar boson channels. 
First, the 14 neutral channels are expressed in the weak eigenbasis as 




), -^\hihi), \hiZi), \z 1 Z 2 ), \h 1 h 2 ), \h 1 Z 2 ), \h 2 Z 1 ), (26) 


where i = 1,2. By taking an appropriate basis transformation, we obtain the following block- 
diagonalized S-wave amplitude matrix: 


Oq — 

where each submatrix is given by 

/ 3Ai 


IGtt 


^4x4 000 

0 ^ 4 x 4 0 0 

0 0 Z3X3 0 


^4x4 — 


2 A 3 -|- A 4 


\ 3 V 2 A, 


V 0 

0 

0 

x 3 / 


2A3 -|- A4 

3 V 2 Af 

3V2XI 

\ 

3A2 

3 V 2 Af 

3 V 2 X\ 


3 V 2 Af 

A3 -|- 2A4 -|- 3 A^ 

3 Af 


CO 

3 A^ 

A3 -|- 2A4 — 

3 Afy 


^4x4 


■^3x3 


( 

X4 

V 2 Xi 

v^A' ) 

A4 

A2 

V 2 Xf 

V 2 X{ 

V 2 X^ 

V 2 Xf 

A3 + A^ 


\V^Xq 

V 2 A{ 

K 

As-Ag^j 


Ai V^(Ag^ -|-iAg) 

Af - iXi A 2 - iX\) 

\V^{Xq — iAg) '/ 2 {Xf + iXlj) A3 -|- A4 j 


(27) 


(28) 


(29) 


(30) 
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Each of the submatrices is obtained in the following basis: 


(^4x4) = 


1 


, z-fZ-, hi hi 

U+UjV + 


V 2 

1 


1 


’ V 2 


I Zf)Zo htyhn 

^ 2^2 + ^ ^ 


OjfoJ2 +002^^1 +-21-^2 + ^ 1 ^ 2 ) 5 o ~ ^1-2:2 + ^2-2^1) J (31) 


1 


^i(^4x4) 


u:{ u}^ — 


z^z-^ 


’ V 2 


^ 2^2 ^ 2^2 


- I —iO{ U)2 — CJo W 


^ + 2 : 4^2 + ^ 1 ^ 2 ) ) ol*^i^‘^2 ~hiZ 2 + h 2 Zi), (32) 


^nzsxs) 


2 2 +^^ 1 ^ 7 ’ ^ 


ZnZn hnhn , 


- 1 + 212:2 - ^ 1^2 + *^ 12^2 + *^ 2 Zi) , 


<{Z3X3) = 


1 

V 2 

1 




.M 2 + M 2 ^ 


2^2 / ! 


— 21^2 + hih 2 + ihi 22 + ih 2 Zi) , 


(33) 


(34) 


where the indices i and j (k and 1) run over 1-4 (1-3). Each of the above bases give the submatrix 
indicated in the parenthesis. We note that at the high energy limit, the hypercharge Y, the isospin 


I and the third component of the isospin I 3 are used to classify the two body scattering states |39], 
namely $a x *1*6 with E = 1 and <I>a x (^^ = iT 2 ^l) with E = 0. In fact, the above bases are 
obtained by finding the two scalar states which belong to the same set of the quantum numbers, 
i.e., the states in the 'I'j^(E 4 x 4 ), 'I'^(E 4 x 4 ) and '^(^^(^ 4 x 4 ) bases respectively belong 

to the state with (E, I, I^) = (0,0,0), (0,1,0), (1,1, —1) and (1,1, —1). 

Next, eight (positive) singly-charged channels are expressed as: 


\ujfzi), \ujfhi), \uj^Z2), la;)^/i 2 ), |wj/ii). 

We obtain the following block diagonalized S-wave amplitude matrix: 

/ El,. n 0 ^ 


(35) 


On - 


1 

IGvr 


>4x4 0 

0 ^3x3 0 

y 0 0 A 3 — A 4 y 

Each of the submatrices and the eigenvalue A 3 — A 4 are obtained in the following basis: 


(36) 


1 


1 


4'f(E4x4) \-iujt 2i +( 2 ^/ 11 ) , —= \ -iuj^Z2+ujth2) 


V 2 

- |—22 — 2 i + uf /i 2 + 0J2 hi) , - |—22 + wj21 — /i 2 + iuj 2 hi) , 


V 2 

1 

2 


1 


1 


1 


(^3x3) =;^ 1*2^1"2i + hi) , ^ |iw^22+a;^h2) , - | 22 + 2i + h2 + cu J hi) 


V 2 


1 


^*^(A 3 - A4) =- |-ia;y22 + iu;^Zi -U]fh2 +uj^hi) , 


(37) 
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where i = 1-4 and j = 1-3. Similar to the discussion for the neutral states, the states in the 
bases 'I'p(y4x4), '^^(Zsxs) and — A4) respectively correspond to the states with {Y, I, 

13 ) = ( 0 , 1 , 1 ), ( 1 , 1 , 0 ) and (1,0,0). 

Finally, three (positive) doubly-charged channels are expressed as: 

1 

We obtain 

We note that the negative charged states are obtained by taking charge conjugation for each positive 
charged channel, which give the same set of eigenvalues of the matrix for the corresponding positive 
charged channels. 

When we impose symmetries in the potential, we obtain the block-diagonalized S-wave ampli¬ 
tude matrix with a smaller size of submatrices. For example, if we assume the CP-invariance (for 
the case of A 5 = Ag = Ay = 0), the maximal size of the submatrices reduces into 3 x 3, or if we 
impose the (softly-broken) Z 2 symmetry (for the case of Ae = A 7 = 0 ), the maximal size of the 
submatrices reduces into 2 x 2 as they have already known in Refs. |p(^ , phf . 

In order to constrain the parameters in the potential, we impose the following condition for 
each eigenvalue of the S-wave amplitude matrix: 




(38) 



|Re(xi)l<^, i = l,...12, (40) 

where ^ is conventionally taken to be 1/2 Q or 1 |Q, and 

Xi = {Eigenvalues(X 4 x 4 ), Eigenvalues(y 4 x 4 ), Eigenvalues(Z 3 x 3 ), A 3 — A 4 }, (41) 

with all Xi being real due to the hermitian nature of the S-wave amplitude matrix. We call the 


bound given by the inequality (^) as the unitarity bound. We note that the analytic expressions 


for the eigenvalues in Eq. (^) are given by solving the fourth and third order equations. In general, 
the solutions of such an equation is given as a too complicated form to explicitly show in this letter. 
Therefore, we do not show the explicit formula, and we numerically calculate the eigenvalues in 
the numerical study given in the next section. 


IV. NUMERICAL STUDIES 

We here discuss the constraint on the parameter space using the unitarity bound. The unitarity 


bound in Eq. (40) sets upper limits on Xj (i = 1,..., 12) given in Eq. (41), which are expressed 
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by combinations of A parameters. Using Eqs. (p^ ) -(|^ ) , this constraint can be translated into the 
bound on physical parameters such as the masses of Higgs bosons and mixing angles. 

It is worthful to mention here that we can obtain an upper limit on the masses of extra Higgs 
bosons even when the SM-like Higgs boson h (= Hi) coupling with the gauge bosons {hVV, 
V = W,Z) slightly deviates from the SM prediction. This can be intuitively understood as follows. 

First, in the Higgs basis defined in Eq. (0), the kinetic terms for the Higgs fields are given by 

^kin = (42) 


where is the covariant derivative. The Higgs-Gauge-Gauge type vertex, i.e., h'lVV, only comes 
from the first term of Eq. (42) at the tree level. In the mass eigenbasis, the ratio of the hVV 
coupling gf^YY in the THDM to that of SM is then expressed by using the rotation matrix R given 


in Eq. (15) as 


„THDM 

^ ^ yhvv _ p 

— SM — 

yhvv 


(43) 


We note that Rn corresponds to sin(/3 — a) in the CP-conserving case as it is seen in Eq. (Al). 
Thus, the non-zero deviation in the hVV couplings from the SM prediction comes from the mixing 
effect of neutral Higgs bosons, i.e., Rn 7 ^ 1. 

Second, when there is no mixing among the neutral Higgs bosons, the mass of h and those of 
the extra Higgs bosons H 2 and are schematically expressed as and M'^ + Xjv'^, respectively. 


as it is seen in Eqs. (1^), (1^) and (1^)- Therefore, in the no mixing case, the upper bound on 
the masses of extra Higgs bosons cannot be obtained, because they can be taken to be as large as 
possible by using the dependence. In other words, we can take the decoupling limit [^] of the 
extra Higgs bosons by taking 

On the other hand, if there is non-zero mixing among the neutral Higgs bosons, i.e., 7 ^ 1, 

dependence appears in which must be kept to be about (125 GeV)^. Therefore, we cannot take 
a too large value of in that case, because we need a large cancellation of the contribution 
to by the AjU^ term which must be excluded by the unitarity bound. Therefore, we can obtain 
an upper limit on the masses of extra Higgs bosons as long as the hVV coupling is deviated from 
the SM prediction. 


In the following, we numerically show the bound on the mass of the second lightest Higgs boson 
denoted as M 2 nd by hxing = 125 GeV and Ky. We assume that the SM-like Higgs boson h is 
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FIG. 1: The upper limit on the mass of the second lightest Higgs boson M 2 nd from the unitarity bound in 
the softly-broken Zi symmetric case. The left and right panels show the Ky and tan /3 dependences of the 
upper limit for fixed values of tan/3 = 1 (black), 3 (blue), 5 (green) and 10 (red) and those of K,y = 0.99 
(black), 0.97 (blue) and 0.95 (red), respectively. The solid (dashed) curve shows the case with ^ = 1/2 (1) 
in Eq. (|^. In both panels, we take (= M 2 nd), and scan the value of M^. 




FIG. 2: The upper limit on the mass of the second lightest Higgs boson M 2 nd from the unitarity bound as 
a function of (= \^') in the case without the GP-violating phases. The black, blue and red curves show 
the case of Ky = 0.99, 0.97 and 0.95, respectively. We take (= M 2 nd) and tan /3 = 1. The 

left (right) panel shows the case of M = {M = m^±/2). 


the lightest of all the Higgs bosons, so that M 2 nd is defined by 

Af 2 nd = hlin(mj|/±, mjy), for the CP-conserving case, 

M 2 nd = hlin(mj|^±, for the CP-violating case. (44) 

We first consider the softly-broken Z 2 symmetric and CP-conserving case, i.e., Ag,? = 0 and 
05 = 0. In this case, we have six free parameters mjj, M^, tan/3 and sin(/3 — a) (= Ky). 
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FIG. 3: Scatter plots for the allowed parameter space from the unitarity bound on the 0cp{= ^6 = 0?) and 
M 2 nd plane in the case of |A 6 |(= IA 7 I) = 0.2 (upper-left), 0.5 (upper-right), 1 (lower-left) and 2 (lower-right). 
For all the plots, we take = rh\ = mjj = M^, 0^=0 and tan /3 = 1, and scan the and sin(/3 — a) 
parameters. The black (circle), blue (square) and red (triangle) dots show the allowed points for the cases 
of Ky = 0.98 ± 0.01, 0.96 ± 0.01 and 0.94 ± 0.01, respectively. 


In Fig. |l], we show the upper limit on M 2 nd as a function of Ky (left panel) and tan/3 (right panel) 
in the case of rnfj± = (= M 2 nd)- In the left (right) panel, each of the curves show 

the cases of tan/3 {ny) to be 1, 3, 5 and 10 (0.99, 0.97 and 0.95). The value of is scanned 


in the both figures. The solid (dashed) curves show the case with ^ = 1/2 (1) in Eq. (|40|) for 
the comparison of the strength of the unitarity constraint with the two cases. Hereafter, we take 
^ = 1/2. From these figures, we can see that the stronger limit on M 2 nd is obtained when 1 — Ky 


is taken to be larger values. Our results for this case are consistent with those given in Ref. |^]. 

Next, we consider the case without the softly-broken Z 2 symmetry but with CP-conservation, 
i.e., ^ 5 , 6,7 = 0. In Fig. we show the upper limit on M 2 nd as a function of \q (= Xj) in the case 
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of tan/3 = 1 and mH± = (= M 2 nd)- The value of is taken to be (left 

panel) and = (M2nd/2)^ (right panel). In both the panels, the black, blue and red curves show 

the cases with Ky = 0.99, 0.97 and 0.95, respectively. Clearly, the stronger bound is obtained in 
the case of = (M 2 nd/ 2 )^ (right panel) than the case of where the appearance 

of kink at around = 0.5, 0.8 and 1 is due to an interchange of eigenvalues which break the 
condition of Eq. (|40|). When is taken to be larger than about 2.2, there is no solution to satisfy 
the unitarity bound. We conhrm that the maximally allowed value of M 2 nd is obtained in the case 
of m2 ~ 

Finally, we consider the most general case with CP-violation. In this case, we have eleven input 
parameters with v = 246 GeV and rrih = 125 GeV as shown in (25). In order to satisfy rufi = 125 
GeV, we scan the rhh parameter for each fixed value of the other parameters. In Fig. |3|, we show 
the allowed parameter regions from the unitarity bound on the 0cp(= ^6 = ^ 7 ) and M 2 nd plane 
in the case of tan/3 = 1, = m\ = mjj = M^ and 65 = 0. The upper-left, upper-right, 

lower-left and lower-right panels respectively show the cases with [Ael = IA 7 I = 0.2, 0.5, 1 and 2. 
In all the panels, the circle (black), square (blue) and triangle (red) points satisfy the unitarity 
bound and Ky = 0.98 ± 0.01, 0.96 ± 0.01 and 0.94 ± 0.01, respectively. The values of and 
sin(/3 — a) are scanned. We find that, in addition to the upper limit on M 2 nd: there is the lower 
limit especially in the case with 9cp 7 ^ 0. The lower bound becomes higher when we take a larger 
value of lAgl (= IA 7 I). The appearance of the lower limit can be understood in the following 
way. If we take a non-zero value of 6 cp, h gives a non-zero value of the off-diagonal mass matrix 
elements and Mig (see Eq. (|T^), which gives non-zero mixings and/or mass splittings among 
the neutral Higgs bosons. On the other hand, we now fix the value of Ky, so that it restricts the 
possible amount of the mixing, and it also requires a non-zero mass splitting among the neutral 
Higgs bosons. Because the SM-like Higgs boson h which we suppose the lightest of all has the mass 
of 125 GeV, the non-zero mass splitting turns out to be the lower limit on M 2 nd- 

We here comment on the constraint from electric dipole moments (EDMs) on the parameter 
space in the THDM. As it has been well known that if a model has an additional source of the 
GP-violation, its magnitude is constrained by EDMs. In Refs. |47, 48], the constraint on parameter 
space from EDMs has been investigated in the softly-broken Z 2 symmetric THDMs. In the THDMs, 
EDMs constrain the allowed region of 9cp depending on the parameter set. Recent study on the 
collider phenomenology of the GP-violating THDM is found in Ref. [49], in which the constraint on 
the complex parameter such as A 5 from EDMs and the unitarity bound turns out to be important. 
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V. CONCLUSIONS 

We have investigated unitarity bounds in the most general two Higgs doublet model without 
a discrete Z 2 symmetry nor CP conservation. We have computed the S-wave amplitudes for two- 
body elastic scatterings of the NG bosons and physical Higgs bosons at high energies for all possible 
initial and final states. By choosing the appropriate bases, the scattering amplitude matrix is given 
to be the block-diagonalized form, and thus the eigenvalues can be easily evaluated numerically. 
We have constrained the parameter space of the model by using the unitarity bound. By fixing 
the mass of the discovered Higgs boson h to be 125 GeV and assuming a non-zero deviation in 
the hVV couplings from the SM values, there is an upper limit on the mass M 2 nd of the second 
lightest Higgs boson. Therefore, by using the precisely measured hVV couplings at future collider 
experiments, we can constrain the allowed region of M 2 nd if a deviation of the hVV coupling from 
the SM limit is found even without its direct discovery. Our results can be useful to constrain 
parameter spaces whenever one evaluates physics quantities in the most general THDM. 

Acknowledgments 

This work was supported, in part, by Grant-in-Aid for Scientihc Research No. 23104006 (SK) 
and Grand H2020-MSCA-RISE-2014 no. 645722 (Non Minimal Higgs) (SK). This work was also 
supported by a JSPS postdoctoral fellowships for research abroad (KY). 

Appendix A: Mass spectrum in specific cases 

We present the mass formulae of the scalar bosons for specific cases of the THDM, i.e., the case 
with CP-conservation but without the Z 2 symmetry, and that with the unbroken Z 2 symmetry 
usually referred as the inert doublet model. 

1. The most general case without CP-violating phases 

When we take the limit of A 5 g 7 —>• 0, the CP symmetry is restored in the Higgs potential, and 
the potential is then described by the ten parameters. The mass of is unchanged and is given 
in Eq. (|ll|). 

For the mass matrix of neutral Higgs bosons, the 3x3 matrix given in Eq. (^) becomes a 
block-diagonalized form, i.e., the 2x2 plus 1x1 form. The first 2x2 part with the basis of 
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corresponds to the mass matrix for the CP-even Higgs states, and the remained 1x1 with the 
basis of h'^ (= A) does the squared mass of the CP-odd Higgs boson. For the CP-even states, we 
obtain the mass eigenstates by introducing the mixing angle /3 — a by 


5/3-a V ^ V (Al) 

^l^—a j y H j 

where h is dehned as the SM-like Higgs boson. 

The squared masses of A is given by 



r^A = Mis- 


(A 2 ) 


Masses for H and h and the mixing angle (3 — a are respectively given by solving Eqs. (^), (p^ 
and ( 1 ^) with the replacement {a,rn\j,m\) —>■ 

m\ = MisCp_a + Mi2s‘p_a + ; (A3) 

= M'iis'p-a + Mi 2 (?p-a ~ 2Aff2'S/3-Q.Cy3_Q,_ , (A 4 ) 

= (A5) 

Consequently, eight input parameters can be chosen as follows 

{mjj, m\, rn\j±, M^, tan/3, sin(/3 —a), A^, A^}, (A 6 ) 


and experimental values of v and ruh- In this parameter choice, the five parameters Ai _4 and 
are derived by Eqs. (16)- (p3|) with the replacement {a,m\,m\). 

The mass formulae in the well-known softly-broken Z 2 symmetric case with the CP-invariance 
are easily obtained by taking Xq = Xj = 0, which agree with the formulae given, e.g., in Ref. |32]. 


2. The inert doublet case 

When we impose an exact Z 2 symmetry in the potential, where the two doublets are transformed 
as <l>i —)• -|-<hi and $2 —^ —*^ 2 , the potential is given for the case with = Ae = Ay = 0 in Eq. (ED- 
In this case, the phase of A 5 can be dropped without loss of generality, so that CP-violation does 
not occur. The VEV of ^2 rnust be assumed to be zero to avoid the spontaneous breakdown of 
the Z 2 symmetry, otherwise the m|<h|<I >2 -|- h.c. terms are generated through the A terms, and the 
inert nature, i.e., without couplings to SM particles, is then lost via the mixing with the active 
Higgs field, namely, the SM-like Higgs boson. 





16 


The squared masses of the scalar bosons are calculated by 

y2 y2 y2 

m\ = AiU^, = 777-2 + -^A 345 , 777,^ = 7772 + ”^("^3 + A 4 — ~ ^"2 + (^7) 
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